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Abstract 

Our presentation centers on the consequences of s-channel unitarity, mani- 
fested by soft re-scatterings of the spectator partons in a high energy diffrac- 



O ti ve process, focusing on the calculations of gap survival probabilities. Our 



emphasis is on recent estimates relevant to exclusive diffractive Higgs pro- 
duction at the LHC. To this end, we critically re-examine the comparison of 
the theoretical estimates of large rapidity gap hard di-jets with the measured 
data, and remark on the difficulties in the interpretation of HERA hard di-jet 
photoproduction. 

V" 1 . 1 Introduction 

> , 

A large rapidity gap (LRG) in an hadronic, photo or DIS induced final state is experimentally defined 
as a large gap in the rj — <fi lego plot devoid of produced hadrons. LRG events were suggested [1-4] 
as a signature for Higgs production due to a virtual W — W fusion subprocess. An analogous pQCD 
process, in which a colorless exchange ("hard Pomeron") replaces the virtual W, has a considerably 
larger discovery potential as it leads also to an exclusive p + H + p final state. Assuming the Higgs 
mass to be in the range of 100 — 150 GeV, the calculated rates for this channel, utilizing proton tagging 
are promissing. Indeed, LRG hard di-jets, produced via the same production mechanism, have been 
observed in the Tevatron [5-17] and HERA [18-29]. The experimental LRG di-jets production rates are 
much smaller than the pQCD (or Regge) estimates. Following Bjorken [3,4], the correcting damping 
■ factor is called "LRG survival probability". 

The present summary aims to review and check calculations of the survival probability as applied 
^ ■ to the HERA-Tevatron data and explore the consequences for diffractive LRG channels at LHC with a 

focus on diffractive Higgs production. 

We distinguish between three configurations of di-jets (for details see Ref. [13-17]): 

1) A LRG separates the di-jets system from the other non diffractive final state particles. On the partonic 
level this is a single diffraction (SD) Pomeron exchange process denoted GJJ. 

2) A LRG separates between the two hard jets. This is a double diffraction (DD) denoted JGJ. 

3) Centrally produced di-jets are separated by a LRG on each side of the system. This is a central 
diffraction (CD) two Pomeron exchange process denoted GJJG. This mechanism also leads to diffractive 
exclusive Higgs production. 

We denote the theoretically calculated rate of a LRG channel by F gap . It was noted by Bjorken 
[3,4] that we have to distinguish between the theoretically calculated rate and the actual measured rate 

fgap 

fgap = (| S | ) • Fg a p. (1) 

The proportionality damping factor [30-33] is the survival probability of a LRG. It is the probability of 
a given LRG not to be filled by debris (partons and/or hadrons). These debris originate from the soft 
re-scattering of the spectator partons resulting in a survival probability denoted | S spec {s) | 2 , and/or from 
the gluon radiation emitted by partons taking part in the hard interaction with a corresponding survival 
probability denoted | Sbrem(Ay) | 2 , 

(| S(s, Ay) | 2 ) = (| S spec (s) | 2 ) • (| S brem (Ay) | 2 >. (2) 



s is the cm. energy square of the colliding particles and Ay is the large rapidity gap. Gluon radiation 
from the interacting partons is strongly suppressed by the Sudakov factor [34]. However, since this 
suppression is included in the perturbative calculation (see 4.3) we can neglect (| Sbrem(Ay) | 2 ) in our 
calculations. In the following we denote ( | S spec | 2 ) = S 2 . It is best defined in impact parameter space 
(see 2.1)). Following Bjorken [3,4], the survival probability is determined as the normalized integrated 
product of two quantities 

2 fd 2 b\M*(s,b) \ 2 P s (s,b) 

jd 2 b\M H (s,b)\ 2 ' K) 

M H (s,b) is the amplitude for the LRG diffractive process (soft or hard) of interest. P s (s,b) is the 
probability that no inelastic soft interaction in the re-scattering eikonal chain results in inelasticity of the 
final state at (s, b). 

The organization of this paper is as follows: In Sec.2 we briefly review the role of s-channel uni- 
tarity in high energy soft scattering and the eikonal model. The GLM model [30-33] and its consequent 
survival probabilities [35-37] are presented in Sec.3, including a generalization to a multi channel re- 
scattering model [38,39]. The KKMR model [40-44] and its survival probabilities is presented in Sec.4. 
A discussion and our conclusions are presented in Sec.5. An added short presentation on Monte Carlo 
calculations of S 2 is given in an Appendix. 



2 Unitarity 

Even though soft high energy scattering has been extensively studied experimentally over the last 50 
years, we do not have, as yet, a satisfactory QCD framework to calculate even the gross features of 
this impressive data base. This is just a reflection of our inability to execute QCD calculations in the 
non-perturbative regime. High energy soft scattering is, thus, commonly described by the Regge-pole 
model [45,46]. The theory, motivated by S matrix approach, was introduced more than 40 years ago and 
was soon after followed by a very rich phenomenology. 

The key ingredient of the Regge pole model is the leading Pomeron, whose linear i-dependent 
trajectory is given by 

a F (t) = a P (0) + apt. (4) 

A knowledge of otjp(t) enables a calculation of a to t, °ei and ^f 2 -, whose forward elastic exponential 
slope is given by 

B el = 2B + loljpln (J^J . (5) 

Donnachie and Landshoff (DL) have vigorously promoted [47, 48] an appealing and very simple Regge 
parametrization for total and forward differential elastic hadron-hadron cross sections in which they offer 
a global fit to all available hadron-hadron and photon-hadron total and elastic cross section data. This 
data, above Pl = WGeV, is excellently fitted with universal parameters. We shall be interested only 
in the DL Pomeron with an intercept ajp(0) = 1 + e, where e = 0.0808, which accounts for the high 
energy growing cross sections. Its fitted [49] slope value is a'p = 0.25 GeV~ 2 . 



2.1 S-channel unitarity 

The simple DL parametrization is bound to violate s-channel unitarity at some energy since a e \ grows 
with energy as s 2e , modulu logarithmic corrections, while a to t grows only as s e . The theoretical problems 
at stake are easily identified in an impact b-space representation. 

The elastic scattering amplitude is normalized so that 

^ = n\f el (s,t)\ 2 , (6) 
a tot = 4:irlmf e i(s,0). (7) 



input 




b 



Fig. 1: A pictorial illustration of a high energy b-space elastic amplitude bounded by unitarity and analyticity/crossing. In the 
illustration we have an input amplitude which violates the eikonal unitarity bound and an output amplitude obtained after a 
unitarization procedure. 

The elastic amplitude in b-space is defined as 

a el (s,b) = i- J ' dcie^ h f el ( S ,t), (8) 
where t = — q 2 . In this representation 

a tot = 2 J d 2 blm[a el {s,b)}, (9) 

a el = J d 2 b\a el (s,b) | 2 , (10) 

CTin = CTtot - Pel- (11) 

As noted, a simple Regge pole with ajp(0) > 1 will eventually violate s-channel unitarity. The 
question is if this is a future problem to be confronted only at far higher energies than presently avail- 
able, or is it a phenomena which can be identified through experimental signatures observed within the 
available high energy data base. It is an easy exercise to check that the DL model [47, 48], with its fitted 
global parameters, will violate the unitarity black bound (see 2.2) at very small b, just above the present 
Tevatron energy. Indeed, CDF reports [50] that a e i(b = 0, y/s = 1800) = 0.96 ± 0.04. A pictorial 
illustration of the above is presented in Fig.l. Note that the energy dependence of the experimental SD 
cross section [13-17] in the ISR-Tevatron energy range is much weaker than the power dependences 
observed for a e i. Diffractive cross sections are not discussed in the DL model. 



2.2 The eikonal model 

The theoretical difficulties, pointed out in the previous subsection, are eliminated once we take into 
account the corrections necessitated by unitarity. The problem is that enforcing unitarity is a model 
dependent procedure. In the following we shall confine ourselves to a Glauber type eikonal model 
[51]. In this approximation, the scattering matrix is diagonal and only repeated elastic re-scatterings are 
summed. Accordingly, we write 

a a (s,b) = i (l - e-^/ 2 ) . (12) 



Since the scattering matrix is diagonal, the unitarity constraint is written as 



2Im[a el {s,b)} = | a el (s,b) | 2 + G in (s,b), (13) 

with 

G in = 1 _ e -n( S ,6)_ (14) 

The eikonal expressions for the soft cross sections of interest are 

a tot = 2(<fb{l-e-^m, (15) 
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d %(l _ e-n^Vy, (16) 
d* 6 (l _ e-n^)), (17) 
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-Q(s,6)/2 



- 2 j (i _ r w) • (18) 

From Eq. (fm i it follows that P s (s, b) = e~ n ( s,b ^ is the probability that the final state of the two initial 
interacting hadrons is elastic, regardless of the eikonal rescattering chain. It is identified, thus, with 
P s (s,b) ofEq.©. 

Following our implicit assumption that, in the high energy limit, hadrons are correct degrees of 
freedom, i.e. they diagonalize the interaction matrix, Eq. (fT2l is a general solution of Eq.(fT3l as long as 
the input opacity SI is arbitrary. In the eikonal model f2 is real and equals the imaginary part of the iterated 
input Born amplitude. The eikonalized amplitude is imaginary. Its analyticity and crossing symmetry 
are easily restored. In a Regge language we substitute, to this end, s Qjp — ► s ap e~ 2* 7rajp . 

In the general case, Eq. (fT3t implies a general bound, | a e [(s, b) |< 2, obtained when G m = 0. 
This is an extreme option in which asymptotically a to t = &ei [52]. This is formally acceptable but not 
very appealing. Assuming that a e \ is imaginary, we obtain that the unitarity bound coincides with the 
black disc bound, | a e i(s,b) |< 1. Accordingly, 

^ < \. (19) 

O'tot ^ 



3 The GLM Model 

The GLM screening correction (SC) model [30-33] is an eikonal model originally conceived so as to 
explain the exceptionally mild energy dependence of soft diffractive cross sections. It utilized the obser- 
vation that s-channel unitarization enforced by the eikonal model operates on a diffractive amplitude in 
a different way than it does on the elastic amplitude. The GLM diffractive damping factor is identical to 
Bjorken's survival probability. 



3.1 The GLM SC model 

In the GLM model, we take a DL type Pomeron exchange amplitude input in which ajp(0) = 1 + A > 
0. The simplicity of the GLM SC model derives from the observation that the eikonal approximation with 
a central Gaussian input, corresponding to an exponential slope of can be summed analytically. This 
is, clearly, an over simplification, but it reproduces the bulk of the data well, i.e. the total and the forward 
elastic cross sections. Accordingly, the eikonal DL type b-space expression for Q(s, b) is: 



Q(s,b) = v(s)T s (s,b), 



(20) 



where, 
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R 2 (s) = ml + Aol P ln{— ), 



and the soft profile is defined 



r s (s,b) 
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(21) 



(22) 



(23) 



irR 2 (s) 

It is defined so as to keep the normalization J d 2 bT s (s,b) = 1. 

One has to distinguish between the eikonal model input and output. The key element is that the 
power A, and v, are input information, not bounded by unitarity, and should not be confused with DL 
effective power e and the corresponding total cross section. Since the DL model reproduces the forward 
elastic amplitude, in the ISR-HERA-Tevatron range, well, we require that the eikonal model output will 
be compatible with the DL results. Obviously, A > e. In a non screened DL type model with a 



Gaussian profile the relation B e \ = 
due to screening. 

With this input we get 

a tot = 2ttR 2 (s) 



hR 2 (s) is exact. In a screened model, like GLM, B e \ > ^R 2 (s) 
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oc -In (s), 



(24) 

(25) 
(26) 



Ei(x) = Jf^ \dt, and C = 0.5773 is the Euler constant. An important consequence of the above is 
that the ratio is a single variable function of v(s). In practice it means that given the experimental 
value of this ratio at a given energy we can obtain an "experimental" value of v which does not depend 
on the adjustment of free parameters. 

The formalism presented above is extended to diffractive channels through the observation, traced 
to Eqs.(|3j and (fT4)h that P s (s,b) = e~ n( - s ' b \ Accordingly, a screened non elastic diffractive cross 
section is obtained by convoluting its b-space amplitude square with the probability P s . 

The above has been utilized [30-33] to calculate the soft integrated single diffraction cross sec- 
tion. To this end, we write, in the triple Regge approximation [53], the double differential cross section 



M 2 d<j sd 
dhPdt 



, where M is the diffracted mass. We, then, transform it to b-space, multiply by P (s, b) and 



integrate. The output 



changes its high energy behaviour from s 2A modulu Zn( 



f^) (which is 
Note also 



dhPdt 

identical to the behaviour of a DL elastic cross section) to the moderate behaviour of ln(— 
a major difference in the diffractive b-space profile which changes from an input central Gaussian to an 
output peripheral distribution peaking at higher b. Consequently, the GLM model is compatible with the 
Pumplin bound [54, 55]. 



<7el(s,b) + (Jdiff(s,b) 1 



&tot(s,b) 



(27) 



3.2 Extension to a multi channel model 

The most serious deficiency of a single channel eikonal model is inherent, as the model considers only 
elastic rescatterings. This is incompatible with the relatively large diffractive cross section observed 
in the ISR-Tevatron energy range. To this we add a specific problematic feature of the GLM model. 
Whereas, atot, o e i and B e i are very well fitted, the reproduction of a s d, in the available ISR-Tevatron 



range, is poorer. A possible remedy to these deficiencies is to replace the one channel with a multi 
channel eikonal model, in which inelastic diffractive intermediate re-scatterings are included as well [38, 
39,56]. However, we have to insure that a multi channel model does improve the diffractive (specifically 
SD) predictions of the GLM model, while maintaining, simultaneously, its excellent reproductions [30- 
33] of the forward elastic amplitude, as well as its appealing results on LRG survival probabilities [35-37] 
to be discussed in 3.3. 

In the simplest approximation we consider diffraction as a single hadronic state. We have, thus, 
two orthogonal wave functions 

(* h | * d ) = 0. (28) 

is the wave function of the incoming hadron, and ^ is the wave function of the outgoing diffractive 
system initiated by the incoming hadron. Denote the interaction operator by T and consider two wave 
functions and ^ which are diagonal with respect to T. The amplitude of the interaction is given by 

Ai tk = {^i^ k | T | = a i)fc 5 ijit 8 k)k >. (29) 

In a 2 x 2 model i, k = 1, 2. The amplitude a^ k satisfies the diagonal unitarity condition (see Eq. (fT3l ) 

2Ima hk (s,b) =\a hk (s,b) \ 2 + G™ k {s,b), (30) 
for which we write the solution 

ai,k(s,b) = i I 1 - e 2 J , (31) 

and 

Gf fc = l- e - n '^ b l (32) 
rij fc(s, b) is the opacity of the (i, k) channel with a wave function Vl/j x 

SU >k = v i)k ( 8 )rf ik ( 8 ,b) (33) 

where 

A 



"u- <k[j- ) ' (34) 



s 



The factorizable radii are given by 



Rl k (s) = 2Rj fl + 2Rl k + Aa'jplni-). (35) 

so 

TfJs, b) is the soft profile of the (i,k) channel. The probability that the final state of two interacting 
hadron states, with quantum numbers i and k, will be elastic regardless of the intermediate rescatterings 
is 

Pf ;k (s,b) = e -^(^) = {1 - a i>fc ( S ,6)} 2 . (36) 
In the above diagonal representation, ^ and ^ can be written as 

V h = a*i + 0V 2 , (37) 

* d = -/Mi + a^ 2 . (38) 
V&i and ^2 are orthogonal. Since | ^ | 2 = 1, we have 



a 2 + p 2 = 1. 



(39) 



The wave function of the final state is 



=| T | V h x V h ) = 

a 2 ai i i{^'i x ^1} + a/3ai, 2 {*i x * 2 + ^2 X ^1} + 

/3 2 a 2j2 {# 2 x tf 2 }. (40) 

We have to consider 4 possible re-scattering processes. However, in the case of a pp (or pp) collision, 
single diffraction at the proton vertex equals single diffraction at the antiproton vertex, i.e., ai j2 = « 2 ,i 
and we end with three channels whose b-space amplitudes are given by 

a e i(s,b) = (V h x ^ h I */) = a 4 ai 5 i + 2a 2 f3 2 a^ 2 + /? 4 a 2j2 , (41) 

a s d(s, b) = {V h x m d I V f ) = af3{a 2 a 1A + (a 2 - f3 2 )a 1)2 + /3 2 a 2 , 2 }, (42) 

a d d(s,b) = (* d x * d I */) = a 2 f3 2 {a 1A - 2a 1>2 + a 2 , 2 }. (43) 

In the numeric calculations one may further neglect the double diffraction channel which is exceedingly 
small in the ISR-Tevatron range. This is obtained by setting a 2)2 = 2ai j2 — a\\. Note that in the limit 
where << 1, we reproduce the single channel model. 

As in the single channel, we simplify the calculation assuming a Gaussian b-space distribution of 
the input opacities soft profiles 

The opacity expressions, just presented, allow us to express the physical observables of interest as func- 
tions of ui ; i, v\£, Riii R12 an( ^ 0' which is a constant of the model. The determination of these 
variables enables us to produce a global fit to the total, elastic and diffractive cross sections as well as the 
elastic forward slope. This has been done in a two channel model, in which add is neglected [38]. The 
main conclusion of this study is that the extension of the GLM model to a multi channel eikonal results 
with a very good overall reproduction of the data. The results maintain the b-space peripherality of the 
diffractive output amplitudes and satisfy the Pumplin bound [54,55]. Note that since different experi- 
mental groups have been using different algorithms to define diffraction, the SD experimental points are 
too scattered to enable a tight theoretical reproduction of the diffractive data, see Fig|2] 



3.3 Survival probabilities of LRG in the GLM model 

The eikonal model simplifies the calculation of the survival probability, Eq.©, associated with the soft 
re-scatterings of the spectator partons. We can, thus, eliminate the nominator and denominator terms in 
I M H (s, b) 1 2 which depend exclusively on s. In the GLM model we assume a Gaussian b-dependence 
for I M H (s, b) \ 2 corresponding to a constant hard radius R h2 . This choice enables an analytic solu- 
tion of Eq.(|3ll- More elaborate choices, such as dipole or multi poles distributions, require a numerical 
evaluation of this equation. 

Define, 

a n(s) = —t 1 - > 1. (45) 



is 



an{s) grows logarithmically with s. As stated, Eq.© can be analytically evaluated with our choice of 
Gaussian profiles and we get 



^ 2 a H (s)j[a H (s),v(s)] 



s z = ^ ' ^ iL v " (46) 




Fig. 2: Integrated SD data and a two channel model fit. 



where j(a, v) denotes the incomplete Euler gamma function 

7(0,2;) = / X z a - x e~ z dz. (47) 



The solution of Eq. d46b . at a given s, depends on the input values of R h2 , R 2 and v(s). In the 
GLM approach, R H<2 is estimated from the excellent HERA data [57-59] on 7 + p — ► J/^S + p. The 
values of v(s) and R 2 (s) are obtained from the experimental pp data. This can be attained from a global 
fit to the soft scattering data [38]. Alternatively, we can obtain v from the ratio and then obtain the 
value of R? from the explicit expressions given in Eqs. (l24l25l26l) . LHC predictions presently depend 
on model calculations with which this information can be obtained. Once we have determined v{s) and 
aff{s), the survival probability is calculated from Eq.d46b. 

In the GLM three channel model we obtain for central hard diffraction of di-jets or Higgs a survival 
probability, 



J d 2 b (a 4 P? Q?/ + 2a 2 /5 2 P? + pPfr fig* 

Scd(s) = ^ 7 : — 5 : i k — — L . (48) 

/ d?b (a 4 n»/ + 2a 2 (3 2 ^ 2 + (3* 



The hard diffractive cross sections in the (i,k) channel are calculated using the multi particle optical 
theorem [53]. They are written in the same form as the soft amplitudes 

tf* = ^(sfT^b), (49) 

where, 



^ k = ^[jt H ■ (50) 



As in the single channel calculation we assume that T^ k (b) is Gaussian, 



r &(&) = ~l^ e (SD 



Note, that the hard radii Rf 2 are constants derived from HERA J/ty photo and DIS production [57-59]. 
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Fig. 3: A contour plot ofS 2 (lC) against u(s) and cih[s). 



As it stands, a three channel calculation is not useful since a^d is very small and the 3'd channel 
introduces additional parameters which can not be constraint by the meager experimental information on 
&dd [13-17]. In a two channel model Eq. d48b reduces to 



A new, unpublished yet, model [60], offers an explicit S 2 calculation for the exclusive NN — > N + 
LRG + 2 J + LRG + N final state, both in one and two channel eikonal models. We shall comment on 
its output in the next subsection. 

3.4 GLM S 2 predictions 

Following are a few general comments on the GLM calculations of S 2 , after which we discuss the 
input/output features of the single and two channel models. Our objective is to present predictions for 



The only available experimental observable with which we can check the theoretical S 2 predic- 
tions is the hard LRG di-jets data obtained in the Tevatron and Hera. A comparison between data and 
our predictions is not immediate as the basic measured observable is f gap and not S 2 . The application of 
the GLM models to a calculation of f gap depends on an external input of a hard diffractive LRG cross 
section which is then corrected by S 2 as presented above. Regardless of this deficiency, the introduction 
of a survival probability is essential so as to understand the huge difference between the pQCD calcu- 
lated Fg a p and its experimental value f gap . A direct test of the GLM predictions calls for a dedicated 
experimental determination of S 2 . The only direct S 2 information from the Tevatron is provided by a 
JGJ ratio measured by DO [5-7] in which g^^^oo) = 2 - 2 ± °- 8 - This is to be compared with a GLM 
ratio of 1.2 — 1.3 ± 0.4 presented below. 

The survival probabilities of the CD, SD and DD channels are not identical. The key difference is 
that each of the above channels has a different hard radius. A measure of the sensitivity of S 2 to changes 
in v and a# is easy to identify in a single channel calculation which is presented in Fig.3. Indeed, 
preliminary CDF GJJG data [17] suggest that f gap measured for this channel is moderately smaller than 
the rate measured for the GJJ channel. 



c2 



jd 2 b (pfrnft 2 - 2(3 2 (p? A n« 1 2 - if 2 iff 2 2 )) 



(52) 




LHC. 



GLM soft profile input is a central Gaussian. This is over simplified, and most models assume 
a power like dipole or multipole b-dependence of V s (s, b) and T H (s, b). Explicit comparisons [60] of 
S 2 obtained with different input profiles shows a diminishing difference between the survival probability 
outputs, provided their effective radii are compatible. 

Regardless of the attractive simplicity of the single channel model, one should add a cautious 
reminder that the single channel model does not reproduce a s d well since its survival probabilities are 
over-estimated. Consequently, we are inclined to suspect that the S 2 values presented in the table below 
are over-estimated as well. 
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As we noted, the soft input can be obtained from either a model fit to the soft scattering data 
or directly from the measured values of ototi^ei an d R r2 . The first method is denoted F1C and the 
second is denoted D1C. Note that having no LHC data, S 2 JD (D1C), at this energy, is calculated on the 
basis of model estimates for the total and elastic cross sections. The constant hard radius R h2 = 7.2 
is deduced from HERA J/ ^ photoproduction forward exponential slope which shows only diminishing 
shrinkage [57,58]. This is a conservative choice which may be changed slightly with the improvement of 
the Tevatron CDF estimates [61] of R h2 . The two sets of results obtained are compatible, even though, 
S 2 (D1C) is consistently lower than S 2 (F1C). The S 2 output presented above depends crucially on the 
quality of the data base from which we obtain the input parameters. The two sets of Tevatron data at 
1800 GeV have a severe 10 — 15% difference resulting in a non trivial ambiguity of the S 2 output. 

The global GLM two channel fit [38] reproduces the soft scattering data (including SD) remark- 
ably well with (3 = 0.464. Its fitted parameters are used for the soft input required for the S 2 calcu- 
lations. Our cross section predictions for LHC are: o to t = 103.8 mb, o e \ = 24.5 mb, a s d = limb 
and B e i = 20.5 GeV~ 2 . The input for the calculation of S 2 requires, in addition to the soft parameters, 
also the values of vf k and Rf 2 - ^ ne nee ded hard radii can be estimated, at present, only for the CD 
channel, where we associate the hard radii R^ x with the hard radius obtained in HERA exclusive J/ ^ 
photoproduction [57,58] and Rf 2 with HERA inclusive J/^> DIS production [59]. Accordingly, we have 
RH 2 = 7.2 GeV' 2 , and R^ 2 = 2.0 GeV' 2 . We do not have experimental input to determine vf k . We 
overcome this difficalty by assuming a Regge-like factorization crf k /af k = constant. Our predictions 
for the CD survival probabilities are: 6.6% at 540 GeV, 5.5% at 1800 GeV and 3.6% at 14000 GeV. 

These results may be compared with a recent, more elaborate, eikonal formulation [60] aiming to 
calculate the survival probability of a final exclusive N + LRG + 2J(orH) + LRG + N state. These 
calculations were done in one and two channel models. The one channel Sq D predicted values are 14.9% 
at 540 GeV, 10.8% at 1800 GeV and 6.0% at 14000 GeV. These values are remarkably similar to the 
GLM one channel output. In the two channel calculations the corresponding predictions are 5.1%, 4.4% 
and 2.7%, which are marginally smaller than the GLM two channel output numbers. 

In our assessment, the two channel calculations provide a more reliable estimate of S 2 since they 
reproduce well the soft scattering forward data. Our estimate for the survival probability associated with 
LHC Higgs production is 2.5% - 4.0%. 

4 The KKMR Model 

The main part of this section (4.1-4.3) was written by V.A. Khoze, A.D. Martin and M. Ryskin (KMR) 
and is presented here without any editing. 



The KKMR model calculation [40-44] of the survival probabilities is conceptually quite similar 
to the GLM model, in as much as unitarization is enforced through an eikonal model whose parame- 
ters provide a good reproduction of the high energy soft scattering data. However, the GLM model is 
confined to a geometrical calculation of S 2 for which we need just the value of R r2 , without any speci- 
fication of the hard dynamics. This value is an external input to the model. The KKMR model contains 
also a detailed pQCD calculation of the hard diffractive proccess, specifically, central diffractive Higgs 
production. Consequently, it can predict a cross section for the channel under investigation. 

4. 1 KKMR model for soft diffraction 

The KMR description [41] of soft diffraction in high energy pp (or pp) collisions embodies 

(i) pion-loop insertions in the bare Pomeron pole, which represent the nearest singularity generated 
by i-channel unitarity, 

(ii) a two-channel eikonal which incorporates the Pomeron cuts generated by elastic and quasi-elastic 
(with N* intermediate states) s -channel unitarity, 

(iii) high-mass diffractive dissociation. 

The KKMR model gives a good description of the data on the total and differential elastic cross 
section throughout the ISR-Tevatron energy interval, see [41]. Surprisingly, KMR found the bare Pomeron 
parameters to be 

A = a(0) - 1 ~ 0.10, a = 0. (53) 

On the other hand it is known that the same data can be described by a simple effective Pomeron pole 
with [47,48,62] 

a f(t) = 1.08 + 0.25 t. (54) 

In this approach the shrinkage of the diffraction cone comes not from the bare pole (a' = 0), but has 
components from the three ingredients, (i)-(iii), of the model. That is, in the ISR-Tevatron energy range 

"o/ ff " = (0.034 + 0.15 + 0.066) GeV~ 2 (55) 

from the 7r-loop, s-channel eikonalisation and diffractive dissociation respectively. Moreover, eikonal 
rescattering suppresses the growth of the cross section and so A ~ 0.10 > A c g ~ 0.08. 

Since the model [41] embodies all the main features of soft diffraction KMR expect it to give 
reliable predictions for the survival probability S 2 of the rapidity gaps against population by secondary 
hadrons from the underlying event, that is hadrons originating from soft rescattering. In particular, KMR 
predict S 2 = 0.10 (0.06) for single diffractive events and S 2 = 0.05 (0.03) for exclusive Higgs boson 
production, pp — > p + H + p, at Tevatron (LHC) energies. 



4.2 Calculation of the exclusive Higgs signal 

The basic mechanism for the exclusive process, pp — > p + H + p, is shown in Fig. |4] The left-hand 
gluon Q is needed to screen the colour flow caused by the active gluons q\ and q2 . Since the dominant 
contribution comes from the region Aq CD <C Q\ <C M#, the amplitude may be calculated using 
perturbative QCD techniques [40, 63] 

M H ~N j ^fgix^x^Ql^fg^x'^Qlf, 2 ), (56) 

where the overall normalisation constant N can be written in terms of the H — > gg decay width [40,64]. 
The probability amplitudes (f g ) to find the appropriate pairs of t-channel gluons (Q, q±) and (Q, 52) are 
given by the skewed unintegrated gluon densities at the hard scale fi, taken to be 0.62 Mh- Since 
the momentum fraction x' transfered through the screening gluon Q is much smaller than that (2) 




Fig. 4: Schematic diagram for central exclusive production, pp — > p + X + p. The presence of Sudakov form factors ensures 
the infrared stability of the Q t integral over the gluon loop. It is also necessary to compute the probability, S 2 , that the rapidity 
gaps survive soft rescattering. 



transferee! through the active gluons (a/ ~ Qtl\fs <C x ~ Muj\fs -C 1), it is possible to express 
fg(x,x',Q 2 ,fi 2 ), to single log accuracy, in terms of the conventional integrated density g{x) [65-68]. 
The / g 's embody a Sudakov suppression factor T, which ensures that the gluon does not radiate in the 
evolution from Q t up to the hard scale n ~ Mh/2, and so preserves the rapidity gaps. 

It is often convenient to use the simplified form [40] 



f g (x, x , Q t , jl 
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(57) 



which holds to 10-20% accuracy. 1 The factor R g accounts for the single log Q 2 skewed effect [67]. It is 
found to be about 1.4 at the Tevatron energy and about 1.2 at the energy of the LHC. 



4.3 The Sudakov factor 

The Sudakov factor T g (Q t , ji) reads [65, 66, 69] 
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with A = kf/(fJL + kt). The square root arises in d57l because the (survival) probability not to emit any 
additional gluons is only relevant to the hard (active) gluon. It is the presence of this Sudakov factor 
which makes the integration in d56l infrared stable, and perturbative QCD applicable 2 . 

It should be emphasized that the presence of the double logarithmic T-factors is a purely classical 
effect, which was first discussed in 1956 by Sudakov in QED. There is strong bremsstrahlung when 
two colour charged gluons 'annihilate' into a heavy neutral object and the probability not to observe 
such a bremsstrahlung is given by the Sudakov form factor 3 . Therefore, any model (with perturbative or 
non-perturbative gluons) must account for the Sudakov suppression when producing exclusively a heavy 
neutral boson via the fusion of two coloured particles. 

'in the actual computations a more precise form, as given by Eq. (26) of [68], was used. 

2 Note also that the Sudakov factor inside t integration induces an additional strong decrease (roughly as M~ 3 [44]) of the 
cross section as the mass M of the centrally produced hard system increases. Therefore, the price to pay for neglecting this 
suppression effect would be to considerably overestimate the central exclusive cross section at large masses. 

3 It is worth mentioning that the H — > gg width and the normalization factor TV in 1561 is an 'inclusive' quantity which 
includes all possible bremsstrahlung processes. To be precise, it is the sum of the H — > gg + ng widths, with n=0,l,2,... . The 
probability of a 'purely exclusive' decay into two gluons is nullified by the same Sudakov suppression. 



More details of the role of the Sudakov suppression can be found in J. Forshaw's review in these 
proceedings [34]. Here KMR would like to recall that the T-factors in [44,70] were calculated to single 
log accuracy. The collinear single logarithms were summed up using the DGLAP equation. To account 
for the 'soft' logarithms (corresponding to the emission of low energy gluons) the one-loop virtual cor- 
rection to the gg — > H vertex was calculated explicitly, and then the scale fi = 0.62 Mh was chosen 
in such a way that eq. d58t reproduces the result of this explicit calculation. It is sufficient to calculate 
just the one-loop correction since it is known that the effect of 'soft' gluon emission exponentiates. Thus 
d58t gives the T-factor to single log accuracy. 

In some sense, the T-factor may be considered as a 'survival' probability not to produce any 
hard gluons during the gg — > H fusion subprocess. However, it is not just a number (i.e. a numerical 
factor) which may be placed in front of the integral (the 'bare amplitude'). Without the T-factors hidden 
in the unintegrated gluon densities f g the integral d56l diverges. From the formal point of view, the 
suppression of the amplitude provided by T-factors is infinitely strong, and without them the integral 
depends crucially on an ad hoc infrared cutoff. 

4.4 Summary of KKMR S 2 predictions 

A compilation of S 2 values obtained in the KKMR model is presented below: 
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A comparison with the corresponding GLM two channel model is possible only for the available GLM 
CD channel, where, the KKMR output is compatible with GLM. KKMR SD and DD output are compat- 
ible with the corresponding GLM single channel numbers. Overall, we consider the two models to be in 
a reasonable agreement. 

A remarkable utilization of the KKMR model is attained when comparing the HERA [18-27] 
and CDF [8-12, 17] di-jets diffractive structure functions derived for the dynamically similar GJJ chan- 
nels. To this end, the comparison is made between the kinematically compatible HERA Fjj(Q 2 = 
75 GeV 2 , (3) and the CDF F?(< E\ >= 75 GeV 2 , 0). The theoretical expectation is that Fjj(/3), as 
measured by the two experiments, should be very similar. As can be seen in Fig|5j the normalizations of 
the two distributions differ by approximately an order of magnitude and for very small j3 < 0.15 there 
is a suggestive change in the CDF distribution shape. This large discrepancy implies a breaking of QCD 
and/or Regge factorization. Reconsidering, it is noted, that HERA DIS data is measured at a high Q 2 
where the partonic interactions induced by the highly virtual photon are point like and, hence, S 2 = 1. 
On the other hand, CDF GJJ measurement is carried out at 1800 GeV and, as we saw, its survival prob- 
ability is rather small. The convolution between the HERA determined GJJ Fjj((3) and the j3 dependent 
survival probabilities, as calculated by KKMR, provides the F®{(3) distribution corrected for the soft 
rescattering of the spectator partons. This is shown in Fig|5]and provides an impressive reproduction of 
the experimental distribution. We were informed [71] that this analysis was successfully redone with an 
updated HI produced structure function distribution. 

The weak element in the above analysis is that it is crucially dependent on the HI determined 
Fjj(P) distribution. ZEUS has constructed a somewhat different structure function. Clearly, a very 
different experimental determination of Fjj((3), such as been recently suggested by Arneodo [72], will 
re-open this analysis for further studies, experimental and theoretical. 
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Fig. 5: The predictions for the diffractive di-jets production at the Tevatron (lower lines), obtained from two alternative sets of 
HERA diffractive parton distributions I and II, compared with the CDF data (shaded area). The upper lines correspond to the 
Tevatron prediction neglecting the survival probability correction. 



4.5 A Comparison between KKMR and GLM 

The approach of GLM and KKMR to the calculation of forward soft scattering in the ISR-Tevatron range 
are basically similar. Both models utilize the eikonal model assuming different input soft profiles which 
have, nevertheless, compatible effective radii. There are, though, a few particular differences between 
the two sets of calculations: 

1) The GLM model, with a Gaussian soft profile, is applicable only in the forward cone (\t\ < 0.3 GeV 2 ), 
where we have most of the data of interest. KKMR use a multipole power behaviour profile which en- 
ables applicability over a, somewhat, wider t range, \t\ < 0.5 GeV 2 . Note that, the GLM output is not 
significantly changed with a multipole power behaviour profile provided its radii are compatible with the 
Gaussian input [60]. 

2) The GLM input Pomeron trajectory is specified by A = 0.12 and a'p = 0.2. These evolve due to 
eikonalisation to an effective output of e = 0.08 and a'jp = 0.25. Note that, A is obtained in GLM as 
a fitted output parameter. In KKMR, the relatively high input A ~ 0.2 is theoretically tuned by a pion 
loop renormalization resulting in an input value of A ~ 0.1. KKMR have a more elaborate treatment of 
ajp(t) than GLM, resulting, nevertheless, with forward cone output predictions similar to GLM. How- 
ever, KKMR accounts for a somewhat wider t range than GLM and reproduces the t dependence of B e \ 
well. Similar results are obtained in a GLM version [39, 56] in which the soft profile is given by a dipole 
distribution. KKMR can predict a few differential properties of S 2 , which are beyond the scope of GLM. 

3) Both models treat the high mass diffraction with the triple Pomeron formalism [53]. In GLM the final 
SD cross section is obtained by a convolution of the input ^fg- with P s (s, b). In KKMR the treatment 
of the SD amplitude is more elaborate, ending, though, with no detailed SD data reconstruction which is 
presented in GLM. 

4) The LHC predictions of the two models for cross sections and slopes are compatible, with the excep- 
tion of add which is neglected in GLM and acquires a significant KKMR predicted value of 9.5 mb. 

GLM is a geometrical model where both the input hard LRG non corrected matrix element squared 
and the soft elastic scattering amplitude, are approximated by central Gaussians in b-space. This property 
enables us to easily calculate the survival probabilities which depend on v, R 2 and R H in a single 
channel input, and on v^, Rf k and Rf k in a two channel input. As we have noted, the GLM model, on 



its own, cannot provide a calculation of F gap and f gap as it needs the hard radii as an external input. The 
KKMR model is more sophisticated. This is attributed to the fact that the hard diffractive LRG process 
is explicitly calculated in pQCD, hence the non corrected F gap and the corrected f gap and Fjj are model 
predictions. As we have just noted, given the hard diffractive matrix element, the actual calculation of the 
diffractive LRG survival probability damping is almost identical to GLM. Keeping this basic observation 
in mind, it is constructive to compare the features of the two models with a special interest on the input 
assumptions and output differences of the two models. 

The main difference between the two models is reflected in the level of complexity of their inputs. 
GLM soft input is obtained from a simple eikonal model for the soft forward scattering, to which we add 
the hard radii which are derived from the HERA data. KKMR calculations of P s are equally simple. 
The calculation of the hard sector matrix elements are, naturally, more cumbersome. Given HERA 
Fjj(Q 2 ,/3), a Tevatron diffractive Fj? in which < Et > and Q 2 are comparable, can be calculated, 
parameter free, without the need to calculate the hard amplitude. But this is a particular case and, in 
general, the KKMR calculation depends on an extended parameter base, such as the the input p.d.f. and 
pQCD cuts. These input parameters are not constrained tightly enough. 

The elaborate structure of the KKMR model provides a rich discovery potential which is reflected 
in the model being able to define and calculate the dependence of S 2 not only on b, but also on other 
variables, notably [3, and experimental cuts such as the recoil proton transverse momentum. GLM de- 
pends on the hard radii external information obtained from HERA data. It lacks the potential richness 
of KKMR. GLM can serve, though, as a standard through which we can compare different unitarized 
models. Given such a model, we can extract effective values for v, R 2 and R h2 and proceed to a simple 
calculation of S 2 . We shall return to this proposed procedure in the final discussion. 

Even though both GLM and KKMR are two channel models, they are dynamically different. GLM 
two channel formulation relates to the diversity of the intermediate soft re-scatterings, i.e. elastic and 
diffractive for which we have different soft amplitudes a^, each of which is convoluted with a different 
probability Pf k which depends on a different interaction radius R 2 k . In the KKMR model the two chan- 
nels relate to two different dynamical options of the hard process. In model A the separation is between 
valence and sea interacting partons. In model B the separation is between small and large dipoles. The 
two models give compatible results. The key point, though, is that the KKMR opacities f^fc, in the 
definition of Pf k , differ in their normalization, but have the same b-dependence. Regardless of this dif- 
ference the output of the GLM and KKMR models is reasonably compatible. The compatibility between 
GLM and KKMR is not surprising since the explicit KKMR calculation of the hard LRG amplitude is 
approximated relatively well by the GLM simple Gaussian. 

Our final conclusion is that the two model output sets are compatible. The richness of the KKMR 
model has a significant discovery potential lacking in GLM. On the other hand, the GLM simplicity 
makes it very suitable as a platform to present different models in a uniform way, which enables a 
transparent comparison. 

5 Discussion 

As we shall see, at the end of this section, there is no significant difference between the values of otot pre- 
dicted by DL and GLM up to the top Cosmic Rays energies. This is, even though, DL is a Regge model 
without unitarity corrections. The explanation for this "paradox" is that the DL amplitude violations of 
s-unitarity are confined, even at super high energies, to small b which does not contribute significantly to 
otot- Note, though, that grows in DL like s e whereas in GLM its growth is continuosly being mod- 
erated with increasing s (see table in 5.3). The DL model predicts that S 2 is identical to unity or very 
close to it in the DL high-i model where a weak IP IP cut is added. The need for survival probabilities so 
as to reproduce the the experimental soft SD cross section values and the hard di-jets rates, is the most 
compelling evidence in support of unitarization at presently available energies. As such, the study of 



high energy soft and hard diffraction serves as a unique probe substantiating the importance of s-channel 
unitarity in the analysis of high energy scattering. 

5.1 S 2 in unitarized models 

Most, but not all, of the unitarized models dealing with LHC S 2 predictions have roughly the same 
S 2 values. This calls for some clarifications. The first part of our discussion centers on the correlated 
investigation of two problems: 

1) How uniform are the output predictions of different unitarization procedures? 

2) How sensitive are the eikonal calculations to the details of the eikonal model they use? 

We start with two non eikonal models which have contradictory predictions. 

The first is a model suggested by Troshin and Tyurin [52]. In this model the single channel 
unitarity constraint (Eq. (fT3b ) is enforced with an asymptotic bound where Gi n = and \a e i\ = 2 i.e. 
asymptotically, a tot = &ei and P s (s,b) = 1. The parameters of the model are set so as to obtain 
a "normal" survival probability monotonically decreasing with energy up to about 2500 GeV where 
it changes its behavior and rises monotonically to its asymptotic limit of 1. Beside the fact that the 
model has a legitimate but non appealing asymptotics, its main deficiency is that it suggests a dramatic 
change in the systematics of S 2 without being able to offer any experimental signature to support this 
claim. Regardless of this criticism, this is a good example of a proper unitarity model whose results are 
profoundly different from the eikonal model predictions. 

Another non eikonal procedure is Goulianos flux renormalization model [17]. This is a phe- 
nomenological model which formally does not enforce unitarity, but rather, a bound of unity on the 
Pomeron flux in diffractive processes. Note that, the Pomeron flux is not uniquely defined so this should 
be regarded as an ad hoc parametrization. Nevertheless, it has scored an impressive success in repro- 
ducing the soft and hard diffractive data in the ISR-Tevatron range. The implied survival probabilities of 
this procedure are compatible with GLM and KKMR. However, the model predicts suppression factors 
for the diffractive channels which are t-independent and, thus, b-independent. The result is that, even 
though the output diffractive cross section is properly reduced relative to its input, there is no change of 
the output profile from its input Gaussian form. Consequently, the Pumplin bound is violated. We are 
informed that Goulianos plans to improve his model by eikonalizing the output of his present model. 

As noted, there are a few eikonal models on the market [73-80], and their predictions are com- 
patible with GLM and KKMR. Reconsidering the procedure of these calculations, their compatibility is 
not surprising once we translate their input to a GLM format. The GLM eikonal S 2 calculation has two 
input sectors in either a single or a two channel version. They are the soft v and R 2 , and the hard radius 
R h2 . Since the soft input is based on a fit of the soft scattering data base, the potential variance in the 
soft parameters is relatively small. The input hard radius is obtained from either the HERA data or a 
theoretical calculation, be it a pQCD diagram or a Regge model. All in all, this is a reasonably stable 
input. In this context, it is interesting to discuss the eikonal model of Block and Halzen [73], where 
the calculated survival probabilities for Higgs production through W-W fusion are seemingly too high, 
5 s (540) = 27%, 5 2 (1800) = 21% and S" 2 (14000) = 13%. Even though, Higgs production is a CD 
process, the above S 2 values are in agreement with the KKMR calculations of S]- )D with a relatively 
high R H = HGeV^ 2 . In a proper S^p calculation, these high S 2 values correspond to an even 
higher R h2 ~ 20 GeV~ 2 , which is far too high as an estimate of the hard radius of WW — ► H. A 
possible interpretation of Block-Halzen results is to associate them with a soft, rather than a hard, LRG 
CD process. This would couple with the non screened interpretation of CD Higgs through the soft CEM 
model [74,75], which predicts very high S 2 values. Since the CEM model is not screened we may, as 
well, assign a survival probability to its output result. This translates into Sq D = S% h Sq EM , providing 
rather reasonable one channel predictions, 5' 2 7D (540) = 18.9% and S£ D (1800) = 7.2%. 

Obviously, each of the eikonal models, quoted above has its own particular presentation and em- 



phasis. They do, however, have compatible results reflecting the observation that their input translates 
into similar values of u, R 2 and R r2 . 

5.2 Compatibility between HERA and the Tevatron di-jets data 

Much attention has been given recently to the compatibility between the Tevatron and HERA DIS GJJ 
data. The starting point made by KKMR and CDF is that rather than depend on a p.d.f. input to calculate 
Fg a p, we may use, the GJJ di-jets diffractive structure function, F®, inferred from HERA DIS data 
[18-27] and associate it with the Fj? derived from the Tevatron GJJ data. As it stands, this procedure 
ignores the role of the survival probability. Consequently, F® obtained from the Tevatron is an order 
of magnitude smaller than the HERA output [8-12, 17,40-44]. This result led to speculations about a 
possible breaking of QCD or Regge factorization or both. Once the Tevatron di-jets diffractive structure 
function is rescaled by the appropriate survival probability, the compatibility between the Tevatron and 
HERA DIS diffractive data is attained. The conclusion of this analysis is that the breaking of factorization 
is attributed to the soft re-scatterings of the the colliding projectiles. Additional hard contribution to the 
factorization breaking due to gluon radiation is suppressed by the Sudakov factor included in the pQCD 
calculation (see 4.3). 

One should note, though, that the HI determination [18-27] of Fj? is not unique. Arneodo [72] 
has suggested a different F® output based on HERA di-jets data which has a different normalization and 
(3 dependences. Should this be verified, there might well be a need to revise the KKMR calculations. 

The evolution of HERA F® from high Q 2 DIS to Q 2 = di-jets photoproduction has raised 
additional concern with regard to the validity of the factorization theorems [28,29]. This is a complicated 
analysis since one has to be careful on two critical elements of the calculations: 

1) The determination of the ratio between direct and resolved exchanged photon (real or virtual). This is 
a crucial element of the theoretical calculation since survival probability is applicable only to the resolved 
photon component. For very high Q 2 data the hard scattering process with the target partons is direct. 
At Q 2 = there is a significant resolved photon contribution. 

2) For di-jets production there is a big difference between the LO and the NLO pQCD calculated cross 
sections [81-83]. Since the HERA analysis compares the pQCD calculation with the di-jets measured 
cross section the normalization and shape of the theoretical input is most crucial in the experimental 
comparison between the high Q 2 and Q 2 = data. 

On the basis of a NLO calculation, Klasen and Kramer [81, 82] conclude that they can reproduce the 
photoproduction data with S 2 = 0.34, applied to the resolved sector. This survival probability is in 
agreement with KKMR and GLM calculations. 

Regardless of the above, preliminary photoproduction GJJ HERA data [28, 29] suggest that both the 
direct and resolved photon sectors are suppressed at Q 2 = 0. A verification of this observation has 
severe consequences for our understanding of the evolution of the diffractive structure function from 
DIS to photoproduction. It does not directly relate, though, to the issue of soft survival probability 
which apply, per definition, only to the resolved photon sector. The suggested effect in the direct photon 
sector should, obviously be subject to a good measure of caution before being substantiated by further 
independent analysis. 

5.3 Diffraction at energies above the LHC 

We end with a table showing the GLM two channel predictions for energies including the LHC, and up 
to the top Cosmic Rays energies. 
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The, somewhat, surprizing observation is that the GLM calculated total cross sections are compatible 
with the DL simple Regge predictions all over the above energy range. This is a reflection of the fact that 
even at exceedingly high energies unitarization reduces the elastic amplitude at small enough b values 
to be relatively insensitive to the calculation of a tot- On the other hand, we see that a e i becomes more 
moderate in its energy dependence and a e i/atot which is 23.6% at the LHC is no more than 26.4% 
at the highest Cosmic Rays energy, \20TeV. The implication of this observation is that the nucleon 
profile becomes darker at a very slow rate and is grey (well below the black disc bound) even at the 
highest energy at which we can hope for a measurment. A check of our results at the Planck scale shows 
o to t = 1010 mb and the profile to be entirely black, i.e., ^eL. = \. is even more moderate in 
its very slow rise with energy. The diminishing rates for soft and hard diffraction at exceedingly high 
energies are a consequence of the monotonic reduction in the values of S 2 with a Planck scale limit of 
S 2 = 0. This picture is bound to have its effect on Cosmic Rays studies. 

Our LHC predictions are compatible with KMR. Note, though, that: i) cr^ LM is rising slowly 
with s gaining 20% from the Tevatron to LHC. KMR has a much faster rise with energy, where, o'^ MR 
is gaining 77% - 92% over the same energy interval, ii) At the LHC B^ LM = 20.5 GeV~ 2 , to be 
compared with a DL slope of 19GeV~ 2 and a KMR slope of 22GeV~ 2 . The GLM 30 TeV cross 
sections are compatible with Block-Halzen. 
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Appendix: Monte Carlo modeling of gap survival 

The following was contributed by Leif Lonnblad and is presented without any editing. 

An alternative approach to gap survival and factorization breaking is to implement multiple inter- 
actions in Monte Carlo event generators. These models are typically based on the eikonalization of the 
partonic cross section in hadronic collisions and can be combined with any hard sub process to describe 
the additional production of hadrons due to secondary partonic scatterings. Some of these programs, 
such as Pythia [84, 85] and HERWIG/Jimmy [86-88], are described in some detail elsewhere in these 
proceedings [89]. Common for all these models is that they include exact kinematics and flavour conser- 
vation, which introduces some non-trivial effects and makes the multiple scatterings process-dependent. 
Also, the predictions of the models are very sensitive to the cutoff used to regularize the partonic cross 
section and to the assumptions made about the distribution of partons in impact parameter space. Never- 
theless, the models are quite successful in describing sensitive final-state observables such as multiplicity 
distributions and jet-pedestal effects [89]. In particular this is true for the model in Pythia which has 
been successfully tuned to Tevatron data 4 by Rick Field [90], the so-called CDF tune A. 

4 Note that the model in PYTHIA has recently been revised [89]. However, the reproduction of Tevatron data is not as good 
for the revised model. 



The Pythia model does not make any prediction for the energy dependence of the total cross 
section - rather this is an input to the model used to obtain the distribution in the number of multiple 
interactions. PYTHIA can, however, make predictions for gap survival probabilities. This was first done 
for Higgs production via W-fusion [2], and amounts to simply counting the fraction of events which 
do not have any additional scatterings besides the W-fusion process. The basic assumption is that any 
additional partonic scattering would involve a colour exchange which would destroy any rapidity gap 
introduced by W-fusion process. Since Pythia produces complete events, these can also be directly 
analyzed with the proper experimental cuts. A similar estimate was obtained for the gaps between jets 
process, both for the Tevatron and HERA case [91]. 

Recently, Pythia was used to estimate gap survival probabilities also for the case of central ex- 
clusive Higgs production [92]. As in the case of gaps between jets, the actual signal process is not 
implemented in Pythia, so direct analysis with proper experimental cuts was not possible. Instead a 
similar hard sub process was used (standard inclusive Higgs production via gluon fusion in this case) 
and the fraction of events without additional secondary partonic scatterings was identified with the gap 
survival probability. Using the CDF tune A the gap survival probability was estimated to be 0.040 for 
the Tevatron and 0.026 for the LHC. This is remarkably close both to the values used in [64] obtained in 
the KKMR model [43], and to the GLM values presented in section l3~4l especiallv the two-channel ones 
obtained in [60]. 
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